Dual-wield NTPases: A novel protein family mined from AlphaFold DB.

AlphaFold protein structure database (AlphaFold DB) archives a vast number of predicted models. We conducted systematic data mining against AlphaFold DB and discovered an uncharacterized P-loop NTPase family. The structure of the protein family was surprisingly novel, showing an atypical topology for P-loop NTPases, noticeable twofold symmetry, and two pairs of independent putative active sites. Our findings show that structural data mining is a powerful approach to identifying undiscovered protein families.

Elastic network model reveals distinct flexibilities of capping proteins bound to CARMIL and twinfilin-tail.

Capping protein (CP) binds to the barbed end of an actin-filament and inhibits its elongation. CARMIL binds CP and dissociates it from the barbed end of the actin-filament. The binding of CARMIL peptide alters the flexibility of CP, which is considered to facilitate the dissociation. Twinfilin also binds to CP through its C-terminal tail. The complex structures of the CP/twinfilin-tail (TW-tail) peptide indicate that the binding sites of CARMIL and TW-tail overlap. However, TW-tail binding does not facilitate the dissociation of CP from the barbed end. We extensively investigated the flexibilities of CP in the CP/TW-tail or CP/CARMIL complexes using an elastic network model and concluded that TW-tail binding does not alter the flexibility of CP. Our extensive analysis also highlighted that the strong contacts of peptides with the two domains of CP, that is, the CP-L and CP-S domains, are key to changing the flexibilities of CP. CARMIL peptides can interact strongly with both of the domains, while TW-tail peptides exclusively interact with the CP-S domain because the binding site of TW-tail on CP relatively shifts to the CP-S domain compared with that of CP/CARMIL. This result supports our hypothesis that the dissociation of CP from the barbed end is regulated by the flexibility of CP.

Rational design of phase separating peptides based on phase separating protein sequence of p53.

Artificial phase-separating (PS) peptides can be used in various applications such as microreactors and drug delivery; however, the design of artificial PS peptides remains a challenge. This can be attributed to the limitation of PS-relevant residues that drive phase separation by interactions of their pairs in short peptides and the difficulty in the design involving interaction with target PS proteins. In this study, we propose a rational method to design artificial PS peptides that satisfy the requirements of liquid droplet formation and co-phase separation with target PS proteins based on the target PS protein sequence. As a proof of concept, we designed five artificial peptides from the model PS protein p53 using this method and confirmed their PS properties using differential interference contrast and fluorescence microscopy. Single-molecule fluorescent tracking demonstrated rapid diffusion of the designed peptides in their droplets compared to that of p53 in p53 droplets. In addition, size-dependent uptake of p53 oligomers was observed in the designed peptide droplets. Large oligomers were excluded from the droplet voids and localized on the droplet surface. The uptake of high-order p53 oligomers into the droplets was enhanced by the elongated linker of the designed peptides. Furthermore, we found that the designed peptide droplets recruited p53 to suppress gel-like aggregate formation. Finally, we discuss aspects that were crucial in the successful design of the artificial PS peptides.

Prediction of chaperonin GroE substrates using small structural patterns of proteins.

Molecular chaperones are indispensable proteins that assist the folding of aggregation-prone proteins into their functional native states, thereby maintaining organized cellular systems. Two of the best-characterized chaperones are the Escherichia coli chaperonins GroEL and GroES (GroE), for which in vivo obligate substrates have been identified by proteome-wide experiments. These substrates comprise various proteins but exhibit remarkable structural features. They include a number of α/ß proteins, particularly those adopting the TIM ß/α barrel fold. This observation led us to speculate that GroE obligate substrates share a structural motif. Based on this hypothesis, we exhaustively compared substrate structures with the MICAN alignment tool, which detects common structural patterns while ignoring the connectivity or orientation of secondary structural elements. We selected four (or five) substructures with hydrophobic indices that were mostly included in substrates and excluded in others, and developed a GroE obligate substrate discriminator. The substructures are structurally similar and superimposable on the 2-layer 2α4ß sandwich, the most popular protein substructure, implying that targeting this structural pattern is a useful strategy for GroE to assist numerous proteins. Seventeen false positives predicted by our methods were experimentally examined using GroE-depleted cells, and 9 proteins were confirmed to be novel GroE obligate substrates. Together, these results demonstrate the utility of our common substructure hypothesis and prediction method.

Structures and mechanisms of actin ATP hydrolysis.

The major cytoskeleton protein actin undergoes cyclic transitions between the monomeric G-form and the filamentous F-form, which drive organelle transport and cell motility. This mechanical work is driven by the ATPase activity at the catalytic site in the F-form. For deeper understanding of the actin cellular functions, the reaction mechanism must be elucidated. Here, we show that a single actin molecule is trapped in the F-form by fragmin domain-1 binding and present their crystal structures in the ATP analog-, ADP-Pi-, and ADP-bound forms, at 1.15-Å resolutions. The G-to-F conformational transition shifts the side chains of Gln137 and His161, which relocate four water molecules including W1 (attacking water) and W2 (helping water) to facilitate the hydrolysis. By applying quantum mechanics/molecular mechanics calculations to the structures, we have revealed a consistent and comprehensive reaction path of ATP hydrolysis by the F-form actin. The reaction path consists of four steps: 1) W1 and W2 rotations; 2) PG-O3B bond cleavage; 3) four concomitant events: W1-PO3- formation, OH- and proton cleavage, nucleophilic attack by the OH- against PG, and the abstracted proton transfer; and 4) proton relocation that stabilizes the ADP-Pi-bound F-form actin. The mechanism explains the slow rate of ATP hydrolysis by actin and the irreversibility of the hydrolysis reaction. While the catalytic strategy of actin ATP hydrolysis is essentially the same as those of motor proteins like myosin, the process after the hydrolysis is distinct and discussed in terms of Pi release, F-form destabilization, and global conformational changes.

Structural Insights into the Regulation of Actin Capping Protein by Twinfilin C-terminal Tail.

Twinfilin is a conserved actin regulator that interacts with actin capping protein (CP) via C terminus residues (TWtail) that exhibits sequence similarity with the CP interaction (CPI) motif of CARMIL. Here we report the crystal structure of TWtail in complex with CP. Our structure showed that although TWtail and CARMIL CPI bind CP to an overlapping surface via their middle regions, they exhibit different CP-binding modes at both termini. Consequently, TWtail and CARMIL CPI restrict the CP in distinct conformations of open and closed forms, respectively. Interestingly, V-1, which targets CP away from the TWtail binding site, also favors the open-form CP. Consistently, TWtail forms a stable ternary complex with CP and V-1, a striking contrast to CARMIL CPI, which rapidly dissociates V-1 from CP. Our results demonstrate that TWtail is a unique CP-binding motif that regulates CP in a manner distinct from CARMIL CPI.

Crystal structure of human V-1 in the apo form.

V-1, also known as myotrophin, is a 13 kDa ankyrin-repeat protein that binds and inhibits the heterodimeric actin capping protein (CP), which is a key regulator of cytoskeletal actin dynamics. The crystal structure of V-1 in complex with CP revealed that V-1 recognizes CP via residues spanning several ankyrin repeats. Here, the crystal structure of human V-1 is reported in the absence of the specific ligand at 2.3 Šresolution. In the asymmetric unit, the crystal contains two V-1 monomers that exhibit nearly identical structures (Cα r.m.s.d. of 0.47 Å). The overall structures of the two apo V-1 chains are also highly similar to that of CP-bound V-1 (Cα r.m.s.d.s of <0.50 Å), indicating that CP does not induce a large conformational change in V-1. Detailed structural comparisons using the computational program All Atom Motion Tree revealed that CP binding can be accomplished by minor side-chain rearrangements of several residues. These findings are consistent with the known biological role of V-1, in which it globally inhibits CP in the cytoplasm.

Protein kinases phosphorylate long disordered regions in intrinsically disordered proteins.

Phosphorylation is a major post-translational modification that plays a central role in signaling pathways. Protein kinases phosphorylate substrates (phosphoproteins) by adding phosphate at Ser/Thr or Tyr residues (phosphosites). A large amount of data identifying and describing phosphosites in phosphoproteins has been reported but the specificity of phosphorylation is not fully resolved. In this report, data of kinase-substrate pairs identified by the Kinase-Interacting Substrate Screening (KISS) method were used to analyze phosphosites in intrinsically disordered regions (IDRs) of intrinsically disordered proteins. We compared phosphorylated and nonphosphorylated IDRs and found that the phosphorylated IDRs were significantly longer than nonphosphorylated IDRs. The phosphorylated IDR is often the longest IDR (71%) in a phosphoprotein when only a single phosphosite exists in the IDR, and when the phosphoprotein has multiple phosphosites in an IDR(s), the phosphosites are primarily localized in a single IDR (78%) and this IDR is usually the longest one (81%). We constructed a stochastic model of phosphorylation to estimate the effect of IDR length. The model that accounted for IDR length produced more realistic results when compared with a model that excluded the IDR length. We propose that the IDR length is a significant determinant for locating kinase phosphorylation sites in phosphoproteins.

All Atom Motion Tree detects side chain-related motions and their coupling with domain motion in proteins.

Structural changes of proteins are closely related with their molecular function. We previously developed a computational tool, Motion Tree (MT), to compare protein structures and describe structural changes using solely the Cα atoms. Here, we have extended MT to incorporate all heavy atoms to analyze side chain-related (SCR) motions. All Atom Motion Tree (AAMT) was applied to 76 proteins that exhibited a simple domain motion identified by MT. AAMT also detected 921 SCR motions. We examined the coupling of domain and SCR motions and classified the structural changes in terms of coupling. The statistical results indicated that it is common for coupled SCR motions to also couple with the domain motion. The classification correlates properties of domain motions and SCR motions. The AAMT results suggest that a large domain motion with a sizable domain boundary is accompanied by SCR motions composed of more than a single residue, which induces further couplings of SCR motions.

Emulsion Droplets Stabilized by Close-Packed Janus Regular Polygonal Particles.

In Pickering-Ramsden emulsions, the packing structure of the colloidal particles at the liquid-liquid (or liquid-gas) interface significantly affects the structure and behavior of the emulsion. Here, using a series of platelike particles with regular polygonal shapes and Janus amphiphilicity, we created emulsion droplets stabilized by close-packed polygonal particles at the interface. The systematic variation of the particle morphology shows that the geometrical features of the regular polygons in (curved) planar packing dominate over the self-assembled structures. The structures are tessellations of triangular, square, and hexagonal particles at the surface for large droplets and regular tetrahedral, cubic, and dodecahedral particle shells of triangular, square, and pentagonal particles for small droplets, respectively. This work creates the possibility of geometrically designing the structure and functionality of emulsions.

Structural changes of homodimers in the PDB.

Protein complexes are involved in various biological phenomena. These complexes are intrinsically flexible, and structural changes are essential to their functions. To perform a large-scale automated analysis of the structural changes of complexes, we combined two original methods. An application, SCPC, compares two structures of protein complexes and decides the match of binding mode. Another application, Motion Tree, identifies rigid-body motions in various sizes and magnitude from the two structural complexes with the same binding mode. This approach was applied to all available homodimers in the Protein Data Bank (PDB). We defined two complex-specific motions: interface motion and subunit-spanning motion. In the former, each subunit of a complex constitutes a rigid body, and the relative movement between subunits occurs at the interface. In the latter, structural parts from distinct subunits constitute a rigid body, providing the relative movement spanning subunits. All structural changes were classified and examined. It was revealed that the complex-specific motions were common in the homodimers, detected in around 40% of families. The dimeric interfaces were likely to be small and flat for interface motion, while large and rugged for subunit-spanning motion. Interface motion was accompanied by a drastic change in contacts at the interface, while the change in the subunit-spanning motion was moderate. These results indicate that the interface properties of homodimers correlated with the type of complex-specific motion. The study demonstrates that the pipeline of SCPC and Motion Tree is useful for the massive analysis of structural change of protein complexes.

Interface property responsible for effective interactions of protean segments: Intrinsically disordered regions that undergo disorder-to-order transitions upon binding.

Proteins that lack a well-defined conformation under native conditions are referred to as intrinsically disordered proteins. When interacting with partner proteins, short regions in disordered proteins can undergo disorder-to-order transitions upon binding; these regions are called protean segments (ProSs). It has been indicated that interactions of ProSs are effective: the number of contacts per residue of ProS interface is large. To reveal the properties of ProS interface that are responsible for the interaction efficiency, we classified the interface into core, rim and support, and analyzed them based on the relative accessible surface area (rASA). Despite the effective interactions, the ProS interface is mainly composed of rim residues, rather than core. The ProS rim is more effective than the rim of heterodimers, because the average rASAs of ProS rim, which is significantly large in the monomeric state, provides a large area to be used for the interactions. The amino acid composition of ProSs correlated well with those of heterodimers in both the core and rim. Therefore, the composition cannot explain why the rASAs of the ProS rim are large in the monomeric state. The balance between a small core and a large rim, and the large solvent exposure of the rim in the monomeric state, are the key to the disorder-to-order transition and the effective interactions of ProSs.

Multiple-Localization and Hub Proteins.

Protein-protein interactions are fundamental for all biological phenomena, and protein-protein interaction networks provide a global view of the interactions. The hub proteins, with many interaction partners, play vital roles in the networks. We investigated the subcellular localizations of proteins in the human network, and found that the ones localized in multiple subcellular compartments, especially the nucleus/cytoplasm proteins (NCP), the cytoplasm/cell membrane proteins (CMP), and the nucleus/cytoplasm/cell membrane proteins (NCMP), tend to be hubs. Examinations of keywords suggested that among NCP, those related to post-translational modifications and transcription functions are the major contributors to the large number of interactions. These types of proteins are characterized by a multi-domain architecture and intrinsic disorder. A survey of the typical hub proteins with prominent numbers of interaction partners in the type revealed that most are either transcription factors or co-regulators involved in signaling pathways. They translocate from the cytoplasm to the nucleus, triggered by the phosphorylation and/or ubiquitination of intrinsically disordered regions. Among CMP and NCMP, the contributors to the numerous interactions are related to either kinase or ubiquitin ligase activity. Many of them reside on the cytoplasmic side of the cell membrane, and act as the upstream regulators of signaling pathways. Overall, these hub proteins function to transfer external signals to the nucleus, through the cell membrane and the cytoplasm. Our analysis suggests that multiple-localization is a crucial concept to characterize groups of hub proteins and their biological functions in cellular information processing.

Comprehensive analysis of motions in molecular dynamics trajectories of the actin capping protein and its inhibitor complexes.

The actin capping protein (CP) binds to actin filaments to block further elongation. The capping activity is inhibited by proteins V-1 and CARMIL interacting with CP via steric and allosteric mechanisms, respectively. The crystal structures of free CP, CP/V-1, and CP/CARMIL complexes suggest that the binding of CARMIL alters the flexibility of CP rather than the overall structure of CP, and this is an allosteric inhibition mechanism. Here, we performed molecular dynamics (MD) simulations of CP in the free form, and in complex with CARMIL or V-1. The resulting trajectories were analyzed exhaustively using Motion Tree, which identifies various rigid-body motions ranging from small local motions to large domain motions. After enumerating all the motions, CP flexibilities with different ligands were characterized by a list of frequencies for 20 dominant rigid-body motions, some of which were not identified in previous studies. The comparative analysis highlights the influence of the binding of the CARMIL peptide to CP flexibility. In free CP and the CP/V-1 complex, domain motions around a large crevice between the N-stalk and the CP-S domain occur frequently. The CARMIL peptide binds the crevice and suppresses the motions effectively. In addition, the binding of the CARMIL peptide enhances and alters local motions around the pocket that participates in V-1 binding. These newly identified motions are likely to suppress the binding of V-1 to CP. The observed changes in CP motion provide insights that describe the mechanism of allosteric regulation by CARMIL through modulating CP flexibility. Proteins 2016; 84:948-956. © 2016 Wiley Periodicals, Inc.

Domain Motion Enhanced (DoME) Model for Efficient Conformational Sampling of Multidomain Proteins.

Large conformational changes of multidomain proteins are difficult to simulate using all-atom molecular dynamics (MD) due to the slow time scale. We show that a simple modification of the structure-based coarse-grained (CG) model enables a stable and efficient MD simulation of those proteins. "Motion Tree", a tree diagram that describes conformational changes between two structures in a protein, provides information on rigid structural units (domains) and the magnitudes of domain motions. In our new CG model, which we call the DoME (domain motion enhanced) model, interdomain interactions are defined as being inversely proportional to the magnitude of the domain motions in the diagram, whereas intradomain interactions are kept constant. We applied the DoME model in combination with the Go model to simulations of adenylate kinase (AdK). The results of the DoME-Go simulation are consistent with an all-atom MD simulation for 10 µs as well as known experimental data. Unlike the conventional Go model, the DoME-Go model yields stable simulation trajectories against temperature changes and conformational transitions are easily sampled despite domain rigidity. Evidently, identification of domains and their interfaces is useful approach for CG modeling of multidomain proteins.

Motion Tree Delineates Hierarchical Structure of Protein Dynamics Observed in Molecular Dynamics Simulation.

Molecular dynamics (MD) simulations of proteins provide important information to understand their functional mechanisms, which are, however, likely to be hidden behind their complicated motions with a wide range of spatial and temporal scales. A straightforward and intuitive analysis of protein dynamics observed in MD simulation trajectories is therefore of growing significance with the large increase in both the simulation time and system size. In this study, we propose a novel description of protein motions based on the hierarchical clustering of fluctuations in the inter-atomic distances calculated from an MD trajectory, which constructs a single tree diagram, named a "Motion Tree", to determine a set of rigid-domain pairs hierarchically along with associated inter-domain fluctuations. The method was first applied to the MD trajectory of substrate-free adenylate kinase to clarify the usefulness of the Motion Tree, which illustrated a clear-cut dynamics picture of the inter-domain motions involving the ATP/AMP lid and the core domain together with the associated amplitudes and correlations. The comparison of two Motion Trees calculated from MD simulations of ligand-free and -bound glutamine binding proteins clarified changes in inherent dynamics upon ligand binding appeared in both large domains and a small loop that stabilized ligand molecule. Another application to a huge protein, a multidrug ATP binding cassette (ABC) transporter, captured significant increases of fluctuations upon binding a drug molecule observed in both large scale inter-subunit motions and a motion localized at a transmembrane helix, which may be a trigger to the subsequent structural change from inward-open to outward-open states to transport the drug molecule. These applications demonstrated the capabilities of Motion Trees to provide an at-a-glance view of various sizes of functional motions inherent in the complicated MD trajectory.

Hierarchical domain-motion analysis of conformational changes in sarcoplasmic reticulum Ca²âº-ATPase.

Sarco(endo)plasmic reticulum Ca(2+)-ATPase transports two Ca(2+) per ATP-hydrolyzed across biological membranes against a large concentration gradient by undergoing large conformational changes. Structural studies with X-ray crystallography revealed functional roles of coupled motions between the cytoplasmic domains and the transmembrane helices in individual reaction steps. Here, we employed "Motion Tree (MT)," a tree diagram that describes a conformational change between two structures, and applied it to representative Ca(2+) -ATPase structures. MT provides information of coupled rigid-body motions of the ATPase in individual reaction steps. Fourteen rigid structural units, "common rigid domains (CRDs)" are identified from seven MTs throughout the whole enzymatic reaction cycle. CRDs likely act as not only the structural units, but also the functional units. Some of the functional importance has been newly revealed by the analysis. Stability of each CRD is examined on the morphing trajectories that cover seven conformational transitions. We confirmed that the large conformational changes are realized by the motions only in the flexible regions that connect CRDs. The Ca(2+) -ATPase efficiently utilizes its intrinsic flexibility and rigidity to response different switches like ligand binding/dissociation or ATP hydrolysis. The analysis detects functional motions without extensive biological knowledge of experts, suggesting its general applicability to domain movements in other membrane proteins to deepen the understanding of protein structure and function.

Hierarchical description and extensive classification of protein structural changes by Motion Tree.

The structures of the same protein, determined under different conditions, provide clues toward understanding the role of structural changes in the protein's function. Structural changes are usually identified as rigid-body motions, which are defined using a particular threshold of rigidity, such as domain motions. However, each protein actually undergoes motions with various size and magnitude ranges. In this study, to describe protein structural changes more comprehensively, we propose a method based on hierarchical clustering. This method enables the illustration of a wide range of protein motions in a single tree diagram, named the "Motion Tree". We applied the method to 432 proteins exhibiting large structural changes and classified their Motion Trees in terms of the characteristic indices of the trees. This classification of the Motion Trees revealed clear relationships to their protein functions. Especially, complex structural changes are significantly correlated with multi-step protein functions.

IDEAL in 2014 illustrates interaction networks composed of intrinsically disordered proteins and their binding partners.

IDEAL (Intrinsically Disordered proteins with Extensive Annotations and Literature, http://www.ideal.force.cs.is.nagoya-u.ac.jp/IDEAL/) is a collection of intrinsically disordered proteins (IDPs) that cannot adopt stable globular structures under physiological conditions. Since its previous publication in 2012, the number of entries in IDEAL has almost tripled (120 to 340). In addition to the increase in quantity, the quality of IDEAL has been significantly improved. The new IDEAL incorporates the interactions of IDPs and their binding partners more explicitly, and illustrates the protein-protein interaction (PPI) networks and the structures of protein complexes. Redundant experimental data are arranged based on the clustering of Protein Data Bank entries, and similar sequences with the same binding mode are grouped. As a result, the new IDEAL presents more concise and informative experimental data. Nuclear magnetic resonance (NMR) disorder is annotated in a systematic manner, by identifying the regions with large deviations among the NMR models. The ordered/disordered and new domain predictions by DICHOT are available, as well as the domain assignments by HMMER. Some examples of the PPI networks and the highly deviated regions derived from NMR models will be described, together with other advances. These enhancements will facilitate deeper understanding of IDPs, in terms of their flexibility, plasticity and promiscuity.

Substrate-shielding and hydrolytic reaction in hydrolases.

In general, transferases undergo large structural changes and sequester substrate molecules, to shield them from water. By contrast, hydrolases exhibit only small structural changes, and expose substrate molecules to water. However, some hydrolases deeply bury their substrates within the proteins. To clarify the relationship between substrate-shielding and enzymatic functions, we investigated 70 representative hydrolase structures, and examined the relative accessible surface areas of their substrates. As compared to the hydrolases employing the single displacement reaction, the hydrolases employing the double displacement reaction bury the substrate within the proteins. The exo hydrolases display significantly more substrate-shielding from water than the endo hydrolases. It suggests that the substrate-shielding is related to the chemical reaction mechanism of the hydrolases and the substrate specificity.